Quantum oscillations of self-dual Abrikosov-Nielsen-Olesen vortices 
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The mass shift induced by one-loop quantum fluctuations on self-dual ANO vortices is computed 
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INTRODUCTION 

In this paper we present some new results on the quan- 
tization of the self dual multi-vortex solutions of the 
Abelian Higgs model. We also take the opportunity to of- 
fer a detailed description of the concepts and techniques 
that allowed us to compute the one-loop quantum cor- 
rection to the mass of self-dual Abrikosov-Nielsen-Olesen 
vortices with one quantum of magnetic flux in the Rapid 
Communication dj. The AHM provides a theoretical 
ground in several fields of physics: it provides shape to 
interesting truncations of the electroweak or grand uni- 
fied theories, it also provides the basis for the various 
phenomenological models for cosmic strings, or it can be 
used as a Ginzburg-Landau theory for superconductivity. 

Interest in this research, developed in the supersym- 
metric framework in was rekindled two years ago. 

Non- vanishing quantum corrections to the mass of N = 2 
supersymmetric vortices were recently reported in papers 
i and@, see also [2^. In the second paper, it was found 
that the central charge of the N = 2 SUSY algebra re- 
ceives a non- vanishing one-loop correction that is exactly 
equal to the one-loop mass shift; thus, one could talk in 
terms of one-loop BPS saturation. This result fits in a 
pattern first conjectured in 6] and then proved in [?J for 
supersymmetric kinks. Another work by the authors of 
the Stony Brook/ Viena group, [8j unveils a similar kind 
of behavior of supersymmetric BPS monopoles in N = 2 
SUSY Yang-Mills theory. In this reference, however, it 
is pointed out that (2+l)-dimensional SUSY vortices do 
not behave exactly in the same way as their (1+1)- and 
(3+l)-dimensional cousins. One-loop corrections in the 
vortex case are in no way related to an anomaly in the 
conformal central charge, contrarily to the quantum cor- 
rections for SUSY kinks and monopoles. 

We shall focus, however, on the purely bosonic Abelian 
Higgs model and rely on the heat kernel/generalized 
zeta function regularization method that we developed 
in reference |9j- Our approach profits from the high- 
temperature expansion of the heat function, which is 
compatible with Dirichlet boundary conditions in purely 



bosonic theories. In contrast, the application of a sim- 
ilar regularization method to the supersymmetric kink 
requires SUSY- friendly boundary conditions, see ^(J- I n 
the kink quantum correction in the </> 4 model is es- 
timated by this method and compared with the correct 
answer obtained from the Dashen-Hasslacher-Neveu for- 
mula, |ll| in order to check the reliability of our ap- 
proach. The relative error found is approximately 0.07%. 
In [13 and [3 we also calculated the quantum mass cor- 
rections for kinks arising in two-component scalar mod- 
els, where second-order small fluctuations are ruled by 
matrix differential operators. Therefore, we were led to 
generalize the zeta function method to the matrix case, 
because the DHN approach, based on a direct compu- 
tation of the spectral density, is not efficient for ma- 
trix differential Schrodinger operators. This step has 
proved to be crucial, opening the possibility of apply- 
ing our method to two-dimensional topological defects in 
the Abelian Higgs model. 

In order to accomplish this task we shall en- 
counter more difficulties than for one-dimensional multi- 
component kinks. As noticed by Vassilevich, the lack 
of analytical expressions for vortex solutions forces us to 
perform a numerical analysis already at the classical level 
to solve the field equations. Also, the high-temperature 
expansion of the heat trace becomes more involved due 
to the jump from one to two spatial dimensions; the 
recurrence relations hold between partial -rather than 
ordinary- derivatives of the high-T expansion coefficients. 
We stress that the evaluation of the Seeley coefficients is 
a very laborious task: fluctuations of the vector, Higgs 
and Goldstone fields are governed by one 4 x 4-matrix 
differential operator, whereas fluctuations of the ghosts 
are determined by one scalar differential operator acting 
on L 2 (M. 2 ). There is, however, one point where the sit- 
uation is more favorable as compared to the kink case: 
the generalized zeta function regularization method pro- 
vides us directly with a finite quantity, without the need 
of infinite renormalizations. This fact is peculiar to even 
spatial dimensions and is probably related to the lack of 
anomalies when fermions are added. As for kinks, we 
shall obtain a simple formula for the one-loop quantum 
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mass correction depending on the Seeley coefficients and 
the number of zero modes. 

One remarkable aspect of our results is that the cor- 
rection found by this method in the bosonic system is 
essentially twice the correction arising in the supersym- 
metric case in Q and [jj. This seems to be in agreement 
with the relationship between the supersymmetric and 
non-supersymmetric one-loop corrections to the masses 
of the sine-Gordon and </> 4 kinks, see 

Q and [3. 

The organization of the paper is as follows: In Section 
§.2 we revise the perturbative sector of the Abelian Higgs 
model in the Feynman-'t Hooft renormalizable gauge and 
set the one-loop mass renormalization conventions. Sec- 
tion §.3 is devoted to studying ANO vortex solutions and 
their fluctuations in a partially analytical, partially nu- 
merical manner. The high-temperature expansion of the 
pertinent heat traces is developed in Section §.4 . Sec- 
tion §.5 explains how quantum oscillations of vortices are 
accounted for in the framework of generalized zeta func- 
tion regularization. In Section §.6 the one-loop vortex 
mass shift formula is applied to cylindrically symmetric 
self-dual vortices. We also briefly discuss how the shift 
depends on the distance between centers of a two-vortex 
solution. Finally, we offer a Summary and Outlook. 



Feynman rules in the R-gauge 

The choice of 4> v = 1 as the ground state causes spon- 
taneous symmetry breaking of the Abelian gauge invari- 
ance. In the Feynman-'t Hooft renormalizable gauge, 

R(A^G) = d^-G 

the particle spectrum involves a vector particle Higgs 
and Goldstone scalar particles <p = 1 + H + iG, and a 
complex ghost \. The Feynman rules are read from the 
action, see Reference [lij : 



S = - 

e 



d 3 3 



-M-g^{d a d a + l)}A v 



+ IdpGdOG - \g 2 + -d^Hd^H - ^-H 2 

+ d^ X *d^x-X*X~Y H{H2 + G2) 

+ A^HG - d^GH) + H(A^- X *x) 



-{H 2 + G 2 ) 2 + -(G 2 + H 2 )A^ 



There are four propagators, plus five third-order and five 
fourth-order vertices shown in the next two Tables: 



THE PLANAR ABELIAN HIGGS MODEL 



TABLE I: Propagators 



The model 



Particle Field Propagator- 



Diagram 



The AHM describes the minimal coupling between an 
[/(l)-gauge field and a scalar field in a phase where the 
gauge symmetry is spontaneously broken. Defining non- 
dimensional space-time variables, x M — > and fields, 
4> — ► v(j) = v(4>i + i<f>2), Afj_ — ► vA^, from the vacuum 
expectation value of the Higgs field v and the J7(l)-gauge 
coupling constant e, the action for the Abelian Higgs 
model in (2+l)-dimcnsions reads: 



-V" + \(D li <t>yD»<f>-u{<i>, 



Higgs H(x) 
Goldstone G(x) 
Ghost x 
Vector Boson A 
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v(k 2 
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- 1 
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ieh 




v(k 2 


- 1 


Fie) 










2 -14 


ie) 



with 



K 2 = is the only classically relevant parameter and 
measures the ratio between the square of the masses of 
the Higgs, M 2 = Xv 2 , and vector particles, m 2 = e 2 v 2 ; 
A is the Higgs field self-coupling. We choose a system 
of units where c = 1, but h has dimensions of length 
x mass. Also, we define the metric tensor as; Q^u — - 
diag(l,-l,-l),A*,f = 0,1,2. 



One-loop renormalization 



Defining 



7(c 2 ) 



d 3 k 



(27r) 3 k 2 - c 2 + is 



and bearing in mind that I(k 2 ) = 1(1) + finite part, the 
one-loop divergences in the planar Abelian Higgs model 
can be organized as follows: 
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TABLE II: Third- and fourth-order vertices 



TABLE III: One- loop counter-terms 



Vertex Weight 



Vertex Weight 



-3ik 2 -^- 

ne 



2 V 





he 



he 



he 



^~^he 



-3i/t 2 ^— 
he 



-3i«; 2 — 




He 



2V_ 

he 



2i Ve^ 



V 

2l Ye 9 



Higgs tadpole 



-2z(k 2 + 1)7(1) + finite part 



• Higgs propagator 



-2i(n 2 + 1)7(1) + finite part 



• Goldstone propagator 

+ ( + 



-2i(n 2 + 1)1(1) + finite part 



• Vector boson propagator 



•vxaAAa/wv 



= 2iJ(l) + finite part 

There are no more one-loop divergent graphs. Therefore, 
in a minimal subtraction scheme, we add the diagrams 
shown in the next Table to cancel the divergences in the 
one-loop graphs. 



Diagram 



Weight 



----- -x — 

/ W\jf\A/\/ 



2i(«f + 1)7(1) 
H(k + 1)7(1) 
H{k + 1)7(1) 
-2i7(l) 



This is tantamount to considering that the counter- 
terms 



C 8 



h( K 2 



V, 



c.t. 
A 



-\l(l)A^ 



j(i)[M a -i] (i) 

(2) 



enter into the Lagrangian. 

All the finite parts are proportional to I(k 2 ) — 7(1) 
and they vanish in the critical point between Type I and 
Type II superconductivity, K = 1, to be considered in 
the sequel. Note that the mass of the elementary particles 
for this critical value of k is taken as subtraction point 
so that the counter-terms exactly cancel the divergence 
due to the Higgs tadpole. Therefore, our renormalization 
criterion is equivalent to the renormalization condition 
stipulated in and []| when k 2 = 1. 



ANO SELF-DUAL VORTICES 

Abrikosov-Nielsen-Olesen vortices are topological de- 
fects satisfying the time-independent field equations: 



diFij = J ; \DiD i( t> = 



(3) 



where Jj — | (<fi*Dj<j) — (Dj<f)* 4>) is the electric current. 
They are static and localized solutions for which the en- 
ergy 

E= /d a x[^-i ! tf + ^W)*A0 + |(0*0-l) a ] (4) 

is finite. Thus, ANO vortices comply with the boundary 
conditions on i.e. when r = \J x\ + x\ tends to oo: 

4>*4>\ slo = 1 , Di<j>\ SL = (d i( j> - iA i( t>)\ SL = , (5) 
i.e., (f>\ s i = e ue , I S Z, and A,\ s i = -iffy^s* . 
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First-order equations 

For the value of the coupling constant k 2 — 1, the 
energy functional can be arranged as follows 

E = J ^ (\D 1 <f>±iD 2 <t>\ 2 + [F 12 ± ±(0*0 -l)] 2 )+i| 5 | 

where g = J d 2 xFi 2 = 2nl is the non-dimensional quan- 
tized magnetic flux. Solutions satisfying the first-order 
differential equations 

D 1( j)±iD 2 (j) = ; Fi 2 ±~(c^-1) = 

or, equivalently, 

(Si^i + A x <h) T {d 2 <t> 2 - A 2 0i) = (6) 
±(5 2 0i + A 2 2 ) + (9 1 2 -A 1 1 ) = O (7) 

F 12 ±i(^ + ^-l) = . (8) 

also solve the second-order equations and are called 
ANO self-dual vortices if they also satisfy the boundary 
conditions (J5J. In what follows, we shall focus on solu- 
tions with positive I: i.e., we shall choose the upper signs 
in the first-order equations. 

Self-dual vortices with cylindrical symmetry 

If 6 = arctan^- is the polar angle, the ansatz 

4>i{xi,x 2 ) = f(r)cosl9 , <t>2(xi,x 2 ) = f(r)sinl9 

ck(r) ot(r) 

Ai(xi,x 2 )=—l sin# , A\(xi,x 2 ) = I cos0 

r r 

plugged into the first-order equations @E1IHJ) leads to: 

^ = T^(/ 2 -l) , f=± l -f(r)[l-a(r)] . (9) 
r dr 21 dr r 



Regular solutions of © with the boundary conditions 
lim f(r) — 1, lim a(r) = 1, zeroes of the Higgs and 

r — >oc r — >oo 

vector fields at the origin, /(0) = 0, a(0) = 0, and integer 
magnetic flux, 

g = — J) dxiAi = —I J) l^dxi ^xidx 2 ] _ ^ 

J r— oo J r— oq ^ 

exist and can be found by a mixture of analytical and 
numerical methods. 

Following the procedure developed in ^t|j we obtain 
numerical solutions for the vortex equations 10. Indeed, 
this approach gives the vortex solution in three differ- 
ent ranges of the radial coordinate. For small values of 
r, a power series is tested in the first-order differential 
equations 10, leading to a recurrence relation between 
the coefficients. Reference ^7) also describes the asymp- 
totic behavior of the solutions. Thus, a numerical scheme 
can be implemented by setting a boundary condition in 
a non-singular point of ©, which is obtained from the 
power series for small values of r. This numerical method 
provides us with the behavior of the vortex solutions for 
intermediate distances by means of an interpolating poly- 
nomial which passes through the numerical data. 

The results are shown in figure 1, where the field pro- 
files a(r) and /(r), the magnetic field B(r) = ^r^r and 
the energy density 

*(r) = \{l-f{r)f+ l ^{l-a{r)ff{r) 

are plotted with respect to r for self-dual ANO vortices 
with I = 1, I = 2, I = 3, and I = 4. A three-dimensional 
view of the energy density in the plane is also shown in 
figure 2 for I = 1, I = 2, I = 3, and / = 4 self-dual 
vortices. Note that the / = 1 vortex shows a different 
pattern as compared with flux tubes of several quanta: 
only in the first case is the energy density maximum at 
the origin (the center). 
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Figure 1. Plots of the field profiles a(r) (a) and f(r) (b), the magnetic field B(r) (c), and the energy density e(r) for vortices 
with 1 = 1 (solid line), I = 2 (broken line), I = 3 (broken-doted line) and 1 = 4 (doted line). 
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Figure 2. 3D graphics of the energy density for 1 = 1,1 = 2,1 = 3 and I = 4 self-dual symmetric ANO vortices. 
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Two-vortex solutions with distinct centers 

To tackle the task of building I = 2 ANO self-dual 
solutions formed by two I = 1 vortices with centers sepa- 
rated by a distance d, we follow the work by Jacobs 

I 



and Rebbi. A variational method is implemented in two 
stages: 

In the first stage, trial functions depending only on a 
single variational parameter w are considered: 



<j> u {z,z*) = $(z,z*) 
A"(z,z*) = u> 



U /«(|* - d/2\) + d/2\) + (1 - co) ^ 



\\A) 



z* - d/2 



aV(\z-d/2\) + 



z*+d/2 

r 



2/ 



a^{\z + d/2\)\ +{l-iu)—a ( - 2 \\z\) 



(10) 
(11) 



Here 

z = Xi+ix 2 
and 



A u (z, z*) = Af(z, z*)+iA%(z, z*) 



(d/2) 



z* 2 - (d/2) 2 



is essentially a phase chosen in such a way that the mag- 
netic flux is equal to Air. and stand 
for the functions / and a associated with self-dual solu- 
tions with cylindrical symmetry -obtained in the previous 



subsection- respectively with vorticity / = 1 and 1 = 2. 
Evoking (|10() and Ijllfl we expect that ui = for the case 
id = and uj = 1 for the case d » 1. Plugging (fTUf) and 
(lllfl into the energy functional, we obtain a expression 
E(u)), which is set to be minimized as a function of u>. 

In the second stage the trial functions are refined by 
adding a deformation such that two requirements are ful- 
filled: 1) the scalar field vanishes at the two centers. 2) 
the gauge-invariant quantities associated with the solu- 
tion are symmetric with respect to the reflection z — > z* . 
The invariant ansatz reads: 



J 



N 



(z,z*) = ^( Z ,z*) + $(z,z*)|z 2 -(d/2) 2 |(cosh|zr 1 ^^/ ! 



(zz*) 



i=0 j=0 



A(z,z*) = A^(z,z*) + 



cosh I 



\ ~> \ -» / (zz*y 



i=0 3=0 



\z*) 




* \ J 




z \3 
z* 



*\3 



N 1 (yy*V 
\ " \ " II \ ZZ ) 

i=0 j=0 



(?)' 




z 
z 



Figure 3. 3D graphics of the energy density for I = 2 self-dual separate vortices with centers at distances d = 1, d = 2, d = 3. 
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These expressions involve K = 3 



_ Q (jV+l)(jV+2) 



variational 



parameters /y, <zf,-, ajj. Finding the minimum of the 
energy functional as a function of these H variables - a 
task for Mathematica- a good approximation to the I = 2 
self-dual solution with a distance d between the two I = 



1 vortex centers is obtained. For our purposes setting 
N = 1 such that H — 9 will suffice. The energy density 
for two- vortex solutions found by this method if K = 9 is 
depicted for d = 1, d = 2 and d = 3 in the above figure. 
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Small fluctuations 



We generically denote the vortex solution fields as 
4> V =V = Vi+^2 , A% = V k ,k = l,2 
Assembling the small fluctuations around the solution 

<j>{x) = ip(x) + ip(x) , A k (x) = V k (x) +a k (x) 



( -d 2 

V V>2 



di 



J 



in a four column £(:?), L 2 -integrable second-order fluctu- 
ations around a given vortex solution are still solutions 
of the first-order equations with the same magnetic flux 
if they belong to the kernel of the Dirac-like operator, 
V£{x) = 0, 



-1p2 

-da + Vx 
di + V 2 



1p2 

-di - V 2 
-da + V t J 



( a i(%) \ 
a 2 {x) 
tpi{x) 



r 



The first component of X>£ gives the deformation of the 
vortex equation (jSJ), whereas the third and fourth com- 
ponents are due to the respective deformation of the co- 
variant holomorphy equations (Q) and ©. The second 
component sets the background gauge 

B(a k ,tp; ip) = d k a k - (^1^2 - ^2^1) 



on the fluctuations. 



The operator H^ 
read: 



V^V and its partner H~ = 



H+ = 



/-A+|^| 2 -2Vi^ 2 

-A+|^| 2 -2V 2 V'2 

-2Vi^2 -2V 2 ^ 2 - A +±(3|i/f + 21414-1) 
\ 2Vi^i 2V2V1 ZV k d k 



2V1V1 
2V 2 ^i 
-2Vkd k 
A+i(3|V| 2 + 2\414-l)/ 



n~ 



•A+IVI' 








A+IV-I 2 





A+i(|7A| 2 + l)+1414 
2V k d k 






-2Vkd k 



- A+±(\i>\ 2 + l) + VkV k ) 



One can check that H + arises in the small deforma- 
tion of the second-order equations in the background 
gauge for k = 1, thus ruling the second-order fluctuations 
around the vortex solutions. In fact, for I = one finds 
that H + = H~ = Ho, where 



Ho 



\ 



A - 






1 

-A- 








A- 







■1 
-A+lJ 



is the second-order fluctuation operator around the vac- 
uum in the Feynman-'t Hooft renormalizable gauge: the 
background gauge in the vacuum sector. Note that the 



fluctuations in this gauge correspond to a massive vector 
particle plus scalar Higgs and Goldstone fields. It will be 
useful in the sequel to write the second-order fluctuation 
operators around I > 1 vortices in the form: 



H — Ho 



}t(x)d k + V ± (x) 



where Q k {x) and V r± (x) are 4x4 functional matrices. 
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HIGH-TEMPERATURE EXPANSION OF HEAT 
TRACES 

Index theorem: moduli space of self-dual vortices 

One easily checks that dimkerl?^ = because the 
spectrum of TL~ is definite positive. Thus, the dimen- 
sion of the moduli space of self-dual vortex solutions with 
magnetic charge I is the index of T>: 

ind T> = dim kcr T> — dim ker 

We follow Weinberg [J^, using the background instead 
of the Coulomb gauge, to briefly determine ind T>. The 
spectra of the operators Tt + and Tt~ only differ in the 
number of eigen-functions belonging to their kernels. For 
topological vortices, we do not expect pathologies due to 
asymmetries between the spectral densities of TL + and 
H~, and thus ind I? = Tre"^ - Tre - ^ - . For a 
case in which these asymmetries are important, see the 
treatment of Chern-Simons-Higgs topological and non- 
topological vortices given in |2fJ, |21| . 

The heat trace of a N x N matrix differential operator 

H = H a + Q k (x)d k + V(x) 

-like the Ti^ operators- is defined as 



Tr e 



-0H _ 



tr 



d 2 x Ku{x, x; (3) 



where Kn(x, y\P) is the N x N matrix kernel of the heat 
equation and tr is the usual matrix trace. Therefore, 
Kn{x, y; P) solves the heat equation 



H)K n (x,y;p) = 



dp 

with initial condition 

K n (x,y;0)=l-5W(x-y) 

Because 

p-/3 s 



(12) 



(13) 



K HQ {x,y;p) = — .I-e-^ 

is the heat kernel for the Klein-Gordon operator 7^o, it 
is convenient to write the heat kernel for TL in the form: 

K n (x,mP)=C H (x,mP)K no (a!,tf,p) (14) 

with C H (x,x;0) = 10. Substituting into JT^J we 
find that Cn(x,y; P) solves the transfer equations: 



The high-temperature expansion 

oo 

c n {x,y-,p) = Y,Cn&y-,n)p n 

n=0 

trades the PDE (|15() by the recurrence relations 



[nl+ (xk - Vk){dkl - -^Qk)]c n (x,y;H) = 

- [AI - Qkdk - V]c n ^{x,y-H) (16) 

among the local coefficients with n > 1, with the initial 
condition cq(x,x;H,) = I. Taking into account that 

_ ( g oo 4 . 

Tre"^ = EE / d 2 x[c n } aa (x,x;H)p n = 



n—0 a—1 

-0 00 



(17) 



(15) 



r n=0 

where we have defined the Seeley coefficients as 
c "('^) = y^ J d 2 x[c n ] aa (x,x;H) 



and that the first local coefficient can be easily computed 

ci(x, x; TL) = — V(x) , 

by applying these formulas to the TC operators we obtain 
in the P = -infinite temperature- limit : 

mdV = J-{ Cl (W+)- Cl (W-)} = 

i.e., the dimension of the self-dual vortex moduli space 
is 21. Physically, this means that there are solutions, if 
K = 1, for any location of the /-vortex centers in the 
plane [24| ; all static configurations of self-dual Z-vortices 
can thus be interpreted as states of neutral equilibrium. 

Seeley coefficients 

Computation of the coefficients of the asymptotic ex- 
pansion is a difficult task; to start with, the order two 
local coefficient reads: 

c 2 (x,x;H) = ~ A V(x) + ^Q k {x)Q k {x)V{x) - 

-\d k Qk{x)V{x) + lQk(x)dkV(x) + ]-V 2 {x) . 

D D 2 

Complexity increases strongly for high-order local coeffi- 
cients . 

The recurrence relation l|16|) allows us to express 
Cn(x,y,H) and its derivatives in terms of all the 
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Ck(x,y;Tt) with k < n and their derivatives. One passes 
from this information to the values of the Seeley coef- 
ficients c n (H) in two steps. First, one must reach the 
subtle y —> x limit. In this analytical manoeuvre the 
partial derivatives of c n (x, y;H) at y = x 



play a prominent role. Note also that: 



[c n } ab (x,x;H)^^C- b (x) 



(Ql ' Q2) Cf (f) = Inn 



d^+ a2 [c n ] ab {x,y;H) 
dx^dx* 2 



In the y — ► x limit the recurrence relation (|16|) becomes 



(k + ai + a 2 + l)( ai > a ^Cf +1 (x) = ( ai+2 ' a ^Cf(x) + (Ql ' Q2+2) Cf (x) 



AT ai a 2 / \ / \ r ar+t^ad ar+ts^ad 

u ( gl - r +l,g 3 -f) (jdb fg\ ^_ u ^2 (ai-r.aa-t+l) 



E E E \ ,. M , ) 0.r\i)A ' '' <Ac'7*.c', 

d=l r=0 t=0 x ' 1 ~ 

_^ iV 02 — 1 OL\ 

oE E E a 2 



Cf(af) 



+ 



(18) 



d=l r=0 t=0 
iV aa cti 



t J dx\dx 2 



EEE 

d=l r=0 t=0 



ati \ I ot% 



Qr+tyad 



r J dx\dx 2 



The initial condition cq(x, x] TL) = I means that all 
the (0'T)Qf 6 (£) vanish except (°>°)C# Q (:r) = 1 for a = 
1,2, ■■■,N. Starting from these conditions one com- 
putes all the ^' 7 )C° 6 (5?) local coefficients by using lfl%l) . 
For instance, in order to obtain (°' ^Cg fc (af) for W + we 
need ^' J ^C^ b (x) for /3, 7 = 0,1,2 as data, which in turn 
can be calculated from ^^Cf(x) for 0, 7 = 0, 1, 2, 3, 4, 
and so forth. Evaluation of (°< )q? h (i?) requires knowl- 
edge of 4032 local coefficients !!!. In general, the rule 
is: knowledge of ^ '°^C^ b (x) amounts to knowledge of 
|(n+l)(n + 2)(4n + 3) ^C% h {x) local coefficients with 
k < n. The second step is much simpler: simple numeri- 
cal integration of 5Za=i ( ' ' ' ) C° a (x) over the plane. 



QUANTUM OSCILLATIONS OF SELF-DUAL 
VORTICES 



Standard lore in the semi-classical quantization of soli- 
tons tells us that the one-loop mass shift comes from the 
Casimir energy plus the contribution of the mass renor- 
malization counter-terms: AMy 



Casimir energy and vortex mass renormalization 
counter-terms 



By expanding the static energy Q of the AHM around 
self-dual vortex solutions one obtains, up to second-order 



in £ in the background gauge: 
,.2 



E + — J d x [djdj - ipi(fi2 + ifatpi] 

Also, the ghosts -arising when the quantization proce- 
dure is performed in the background gauge- contribute 
negatively to the energy: 



Thus, the vortex Casimir energy is the sum of the Casimir 
energies of the bosonic 01, a%, <pi, tp 2 fluctuations around 
the vortex minus the Casimir energy of the fermionic 
fluctuation %; the ordinary -non-matrix- Schrodinger op- 
erator ruling the ghost fluctuation around the vortex is: 



ip 2 is a pure gauge oscillation but its contribution is killed 
by the negative ghost contribution. The same applies 
for the vacuum Casimir energy: the Goldstone boson 
Casimir energy is canceled by the ghost Casimir energy, 
the trace of the square root of Hq = — A +1. In sum, 
the vortex Casimir energy measured with respect to the 
vacuum Casimir energy is given by the formal formula: 



AMy 



hm 



STr* (H + ) 2 -STr(Ho) 



9 



STr* (H + ) 2 = Tr* (H + ) 2 - Tr (H G ) 2 

Sir (Wo) 1 = Tr(^o)^ -Tr(^) 

The star means that the 21 zero eigenvalues of TL + must 
be subtracted because zero modes only enter at two-loop 
order. 

In the minimal subtraction renormalization scheme, fi- 
nite renormalizations are adjusted in such a way that the 
critical point k 2 = 1 is reached at first-order in the loop 
expansion. Therefore, and J5J tell us that the contri- 
bution of the mass renormalization counter-terms to the 
vortex mass is: 



and 



AMI 



£(V, v k ) 



AM: 



AM' 



fim/(l)E(^Vfc) 



1, 



dx 2 [{i-i^n-^VkVk} 



and the divergent integral 1(1) can be written in the form 
1(1) = 1 f 1 



(2tt) 2 ^Jk-k + 1 

after applying the residue theorem to integration in the 
complex fco-plane. 



Zeta function regularization of Casimir energies and 
self-energy graphs 

We regularize both infinite quantities AM G and AMy 
by means of generalized zeta functions. From the spectral 
resolution of a Fredholm operator Ti 

one defines the generalized zeta function as the series 

<"<•>- • 

n ™ 

which is a meromorphic function of the complex variable 
s 1^2, |23. We can then hope that, despite their contin- 
uous spectra, our operators fit in this scheme, and write: 

AM G (s) = y(^) {((hHs)-(hc( S )) + 



+ (c*«W-Cwo(*))} 



am*(s) = CH ( S m,v k ) 



where 



Ch (s) 



m 2 L 2 T(s - 1) 

r(s) 



and fi is a parameter of inverse length dimensions. Note 
that 

AMy = lim i AM^(s) , AMy = lim AM^(s) 



7(1) = lim 



1 



2m 2 L 2 



Cn (s) 



on a square of area L 2 . 

Together with the high-temperature expansion, the 
Mcllin transform of the heat trace 



Cn(s) 



1 



r(s) Jo 



s-1 rrs-fm 



d/3/? s_1 Tre 



shows that 



1 f 1 1 

Ch(s) = =pr / df3p s+n ~ 2 c n (H)e^ + —B n (s) 

is the sum of meromorphic and entire —B-n (s)- functions 
of s. Neglecting the entire parts and keeping a finite 
number of terms, No, in the asymptotic series for Ch( s )i 
we find the following approximations for the generalized 
zeta functions concerning the differential operators Ti + 
and Ti G relevant to our problem: 



No 



CnAs)-Cn ( S ) * E c «(^ + ) 



n=l 
No 



7[s + n — 1, 1] 
47rr(s) 



is the incomplete 



7 [.s + n - 1,1] = djSjS- 
gamma function, with a very well known meromorphic 
structure. 

Regarding one-dimensional kinks, see 0, 01 1 Eli 
the contributions of Co(H + ) and co(H G ) to Cw+( s ) an d 
Oh g (s) are respectively canceled by 0h o ( s ) an d Ch g ( s ); 
i.e., renormalization of zero point vacuum energies takes 
care of the c (H + ) and Co(TL G ) contributions to the vor- 
tex Casimir energy. Note, however, that, in contrast to 
the (l+l)-dimensional case, the value s = —\ for which 
we shall obtain the Casimir energy is not a pole. To 
compute the vortex Casimir energy one can first take the 
s = — h limit and then subtract the vacuum Casimir en- 
ergy regularized by this procedure; a finite answer for 
the kink Casimir energy is only reached if one first sub- 
tracts the vacuum Casimir energy of the one-dimensional 
system. 



One-loop mass shift formula 

Writing as c„ = c n (H + ) — c n (TL G ) the difference be- 
tween the Seeley coefficients of Tt + and TL G for vorticity 
I, we check that the contribution of the first coefficient 
to the Casimir energy 



AM^ c (s)^[E- 



Cl 



7[M/2] 
47rr(s) 
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is finite at the s 



■h limit 



AMf c (-l/2) 



7[-l/2,l] 

r(i/2) 



and exactly cancels the contribution of the mass renor- 
malization counter-terms -also finite for s = i— : 



AM*(s) 
AM*(l/2) 



_.E(^V fc )-^j- 



47T 



r(i/2) 



Our choice of a minimal subtraction scheme not only ar- 
ranges finite renormalizations in such a way that self- 
duality holds for k = 1 at the one-loop order, but also 
fits in with the criterion that the mass renormalization 
counter-terms must kill the contribution to the Casimir 
energy of the first Seeley coefficients for the heat trace 
expansions of the operators H + , H G , Ho, Hq . The same 
cancellation happens for kinks only if the mode number 
cut-off regularization procedure, see Q and j2{|, is 
applied. 

Subtracting the contribution of the 21 zero modes, 



AM V = 



Km 



lim 



21 

W)Jo 



dfi(3 s 



N 



j[s + 71—1, 1] 

4tiT(s) 



we finally obtain the following formula for the vortex 
mass shift: 



AAfv 



Hm 



1 ^ r 3 n 



n=2 



21 

7H 



(19) 



ever, depend on successive derivatives of the solution. 
This dependence can increase the error in the estimation 
of these local coefficients because we handle an interpo- 
lating polynomial as the numerically generated solution, 
and the successive derivations with respect to r of such 
a polynomial introduces inaccuracies. Indeed this oper- 
ation is plugged into the algorithm that generates the 
local coefficients in order to speed up this process. It is 
thus of crucial importance to use the first-order differen- 
tial equations JHJ in order to eliminate the derivatives of 
the solution and write the local coefficients as expressions 
depending only on the fields. We find: 



dtpi 
dxi 
chh 
dx 2 
dip 2 
dx\ 
dip 2 
dx 2 
dVi 
dxi 
dVx 
dx 2 
dV 2 
dxi 

dx 2 



r 

lf(r) 
r 

U(r) 
r 

Mr) 



cos 6 cos 19(1 — a(r)) + sin ff sin Iff] 
sin 9 cos 10(1 — a(r)) — cos 6 sin Iff] 
cos 9 sin 19(1 — a(r)) — sin 9 cos Iff] 



r 

sin 9 cos t 
-lcos29 
-lcos29 



s\nffsml9(l — a(r)) + cos 9 cos Iff] 



2lf(r)a(r) + l } 
r 2 

^-ico S ^(/»-l) 

r A 2 

2lf(r)a(r) l., 2/ x 



ONE-LOOP MASS SHIFTS 
Local coefficients for cylindrically symmetric vortices 

We shall apply these formulae to cylindrically sym- 
metric vortices. The heat kernel local coefficients, how- 



for self-dual ANO vortices with generic (positive) vortic- 
ity I. 

The recurrence formula now gives the local coefficients 
of the asymptotic expansion in terms of f(r) and a(r), 
e-g-, 



J 



tv[c 1 ](x,x;H + ) 
tv[c 2 ](x,x;H + ) 

tr[c 3 ](x,x;H + ) 



5[1-/V)]- V 2 «V) 



{37r 4 + 4Z 4 a 4 (r) + 8(7/V - 8r 4 )/ 2 (r) + 27r 4 / 4 (r) 



V-8r 4 )/ 2 ( 

lr 2 a(r)[-l + (1 + 13Z)/ 2 (r)] + 8Z 2 a 2 (r)(-2 - 3r 2 + 9r 2 / 2 (r))} 



1 



{-4i 6 a 6 (r) - 4i 3 r 2 a 3 (r)[14 + (-132 + 1670/ 00] + 4Z 4 o 4 (r)(20 + 9r 2 + 32r 2 / 2 (r)) - 2Zr 2 a(r)[-4(16 + 9r 2 )+ 



120r 6 

+ (64 + 961 - 4721 2 + 344Z 3 + 88Z 2 + 243Zr 2 )/ 2 (r) + (-52 + 109/)r 2 / 4 (r)] + Z 2 a 2 (r)[-256 - 144r 2 - 117r 4 + 



+2r 2 



■ 548Z + 516i 2 + 183r 2 )/ 2 (r) + 99r 4 / 4 (r)] + r 2 [r 2 (-16 + 151r 2 ) + (-320Z 3 + 16(K 4 + 32r 2 + 48lr 2 



-321r 4 + 8Z 2 (20 + 39r 2 ))/ 2 (r) + r 2 (-16 - 48Z + 44Z 2 + 199r 2 )/ 4 (r) - 29r 4 / 6 (r)]} 
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We have explicitly given only the first three local co- 
efficients of the heat kernel expansion for Tt + because 
the complexity of the expressions increases with n enor- 
mously. Additionally 

Cl (x,x;H G ) = l~f(r) 
_1 { [4/ 2 + 5r 2 - 8l 2 a(r) + 4l 2 a 2 (r)\f 2 {r)+ 



c 3 (x,x;H G ) 

+16« 2 (l + r 2 ) 
+16Z 2 (1 - 6Z + 6Z 2 
+16/ 4 a 4 (r)]/ 2 (r) 



g r 2 



1 



-3r 2 -2r 2 / 4 (r)} 
{lOr 4 - [-32/ 3 + 16« 4 +f 



60r 4 

8Z(— 12Z 2 + 8/ 3 + r 2 + 41(1 + r 2 ))a(r) + 
r 2 )o?(r) + 32(1 - 2l)l 3 a 3 (r) + 
r 2 [8l + 1G1 2 + 17r 2 + 16/ 2 Q 2 (r) - 
-8Z(l+4Z)a(r)]/ 4 (r)-4r 4 / 6 (r)} 



TABLE V: Seeley Coefficients for 1 = 3,4. 





I = 3 


1 = 4 


n 


Cn(H+) C n {H G ) 


c n {H+) c n (H G ) 


2 
3 
4 
5 
6 


90.20440 11.51035 
36.68235 2.60898 
11.69979 0.46721 
2.86756 0.067279 
0.566227 0.0079269 


118.67540 16.46895 
46.01141 4.00762 
14.64761 0.77193 
3.58906 0.11747 
0.667202 0.01620 



an enormous number of local coefficients. Nevertheless, 
the choice No = 6 is acceptable. The behavior of the 
asymptotic series in (fB^l is given in Table VI: 



are the first three local coefficients for the heat kernel 
expansion for the ghost operator TL G . 

Plugging these expressions into the partially analytical 
partially numerical solution for /(r) and a(r), it is possi- 
ble to compute the local coefficients and integrate them 
numerically over the whole plane. 



Mass shift for vorticities 1 = 1, 1 = 2, 1 = 3, 1 = 4 

Finally, the one-loop quantum correction of the vortex 
solution with vorticity I is given by formula (|19f) 



AM V = - 



hm 



1 N 
r, /it — ' 



■ill 71 ~ o> x 

n—2 



21 



Using the Mathematica environment in a modest PC we 
have obtained the coefficients shown in Tables IV and V, 



TABLE IV: Seeley Coefficients for 1 = 1,2. 





I = 1 


I = 2 


n 


Cn(H + ) C n {H G ) 


Cn{H + ) C n (H G ) 


2 
3 
4 
5 
6 


30.36316 2.60773 
12.94926 0.31851 
4.22814 0.022887 
1.05116 0.0011928 
0.20094 0.00008803 


61.06679 6.81760 
25.61572 1.34209 
8.21053 0.20481 
2.02107 0.023714 
0.40233 0.002212 



We remark that formula (|19fl depends on the number 
No chosen to cut the asymptotic heat kernel expansions. 
We have no means of determining the optimum value 
for Nq, but in practice we can only cope with a small 
No value; a big No would require the computation of 



TABLE VI: Convergence of the asymptotic series in units of 
hm. 



No 


AMv(No) 


AMy(iVo) 


AMv(No) AMv(N ) 




I = 1 


/ = 2 


I = 3 


I = 4 


2 


-1.02951 


-2.03787 


-3.01187 


-3.97025 


3 


-1.08323 


-2.14111 


-3.15680 


-4.14891 


4 


-1.09270 


-2.15913 


-3.18208 


-4.18014 


5 


-1.09427 


-2.16212 


-3.18628 


-4.18534 


6 


-1.09449 


-2.16257 


-3.18690 


-4.18606 



The convergence up to the sixth order in the asymp- 
totic expansion is very good. In the case of A(0)| kinks 
we found agreement between the result obtained by this 
method and the exact result up to the fourth decimal 
figure, see 0, by choosing N = 10. 

There are reasons to expect this behavior on general 
analytical grounds. Truncation of the asymptotic expan- 
sion of the heat function at order No produces an error of 
order (3 N ° , which in turn leads to an error proportional 
to 7[Vq — 5,1] — N 31 > for Nq large, in the computation 



of the Ch+(— \) ze t a function, see |22j Section 1.10. In 
fact, the rate of convergence is improved in our problem 
by the the smallness of the c„ coefficients, see Tables IV 
and V, for large n. This smallness is due to the fact that, 
when n increases, higher and higher powers of partial 
derivatives of the field profiles of increasing order enter 
in the computation of c n . The vortex solutions, however, 
are as regular and smooth as allowed by the topology. 
Therefore, the admitted error by cutting the mass shift 
formula at No — 6 is especially small for low vorticities. 

In Table VII we give the one-loop quantum corrections 
for the vortex solutions up to I = 4, whereas we plot the 
correction in the figure as a function of the magnetic 
flux. The broken line (linear function) represents the hy- 
pothetical situation in which each magnetic flux quantum 
would contribute with the same correction. Hence, this 
is almost - within the error margin- the situation that we 
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TABLE VII: One-Loop Quantum Mass Correction to the vor- 
tex with vorticity I = 1,2, 3, 4. 



/ 


AM V /Urn 


AM 


1 2 3 4 


1 


-1.09449 


-l 




2 


-2.16257 


-2 




3 


-3.18690 


-3 




4 


-4.18606 


-1 





TABLE IX: Convergence of the asymptotic series. 



iVo 


AM v {N )/hm 


AM v (N )/hm AM v (N )/hm 




d= 1 


d = 2 d = 3 


2 


-2.03770 


-1.99798 -1.98848 


3 


-2.14095 


-2.09695 -2.08672 



SUMMARY AND OUTLOOK 



have found. 

These results, however, do not allow us to answer the 
question of whether or not the classical degeneracy with 
respect to the vortex centers observed at the classical 
level also holds at one-loop order. The figure in Table VII 
seems to suggest that the mass shift of I well separated 
vortices is equal -modulo errors- to I times the mass shift 
of a single vortex, but we do not know in what direction 
the errors run. 



Mass shift for solutions with two separate vortices 

We now offer two Tables, VIII and IX, where Seeley 
coefficients and the quantum corrections are given for 
two- vortex solutions with intermediate separations d = 1 , 
d = 2, and d = 3 between superimposed vortices, w = 
in H10 |I - H11(I . and well separated vortices, ui = 1 in 110(1 - 
iJUJ. The coefficients of the asymptotic expansion are 
computed only up to third order because much more 
computation time is required. Also, we stress that in 
this situation, with no cylindrical symmetry, we expect 
not so good results because there are two more important 
sources of errors: first, the variational solutions with two 
separate vortices are far less exact than the solution with 
/ = 2 and cylindrical symmetry. Second, even though 
another numerical method would be used in the search 
of vortex solution we would run in difficulties; there is 
no way to avoid the use of partial derivatives in the cal- 
culation of the coefficients because the vortex equations 
alone are not enough. 



TABLE VIII: Seeley Coefficients for d = 1, 2, 3. 





d= 1 


d = 2 


d = 3 


n 


Cn(H + ) c n (H a ) 


c„(H + ) c n {H G ) 


Cn(H + ) C n {H G ) 


2 
3 


61.0518 6.81277 
25.6137 1.33822 


58.3359 6.46609 
24.5050 1.23466 


57.3420 6.03872 
24.1187 1.02031 



The one-loop mass shifts of superimposed vortices with 
low magnetic fluxes are: 

(— - 1.09427ft) + o(ft 2 ) 
1.08106ft) + o(ft 2 ) 
1.06230ft) +o(ft 2 ) 
1.04651ft) +o(ft 2 ). 

Much less precise results are also provided for two- 
vortices with separate centers. This is to be compared 
with the supersymmetric result: 



My 


1 


= m ^ 


TTV 

e 


My 


2 


= 2m 


/ 7TV 

\ e 


My= 


3 


= 3m 


/ TTV 






V e 


My= 


4 


= 4m 


("irv 






V e 



M l y 



\l\m - 0.5000ft) +o(ft 2 ) 



see 0] and 0. We notice that the one- loop correction 
due to bosonic fluctuations of self-dual vortices is al- 
most twice the correction arising in the supersymmet- 
ric system coming only from mass renormalization coun- 
terterms when proper SUSY-preserving boundary condi- 
tions are imposed. The same proportion holds between 
one-loop corrections to sine-Gordon and (f) A kink masses 
in the non-supersymmetric and supersymmetric frame- 
works, see 0| and [To] ]. 

It seems plausible that a similar method can success- 
fully be applied to compute the one-loop mass shift for 
self-dual Chern-Simons-Higgs vortices, see H^-|2^. A 
Hamiltonian formalism in the topological sectors of the 
first-order CSH Lagrangian system should be first de- 
veloped. More ambitious, generalized zeta functions of 
12 x 12 matrix PDE operators in three variables are 
essential in computing the one-loop mass shift to BPS 
monopoles. Thus, our procedure opens a door to calcu- 
late quantum corrections to BPS monopole masses in a 
Af = bosonic setting to be contrasted with the N = 2 
and A/" = 4 supersymmetric results of and [2(| . 
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